Ga+, In"*" and T1+ impurities in alkali halide crystals: distortion trends. 
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A computational study of the doping of alkali halide crystals (AX: A = Na, K; X = CI, Br) by 
ns^ cations (Ga^, In""" and Tl""") is presented. Active clusters of increasing size (from 33 to 177 
ions) are considered in order to deal with the large scale distortions induced by the substitutional 
impurities. Those clusters are embedded in accurate quantum environments representing the sur- 
rounding crystalline lattice. The convergence of the distortion results with the size of the active 
cluster is analyced for some selected impurity systems. The most important conclusion from this 
study is that distortions along the (100) and (110) crystallographic directions are not independent. 
Once a reliable cluster model is found, distortion trends as a function of impurity, alkali cation and 
halide anion are identified and discussed. These trends may be useful when analycing other cation 
impurities in similar host lattices. 

PACS numbers; 61.72.Bb; 61.82.Ms 
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I. INTRODUCTION 

Luminescent materials are used in a wide variety 
of technological applications. □ An important proportion 
of these materials is obtained by doping a pure ionic 
crystal with substitutional impurities which have desir- 
able absorption-emission properties. The introduction 
of these substitutional impurities induces a distortion of 
the host lattice in a local region around the defect, due 
mainly to the different sizes of the impurity and host ions. 
That distortion, which is different for each combination 
of impurity ion and host lattice, determines the specific 
lattice potential felt by the impurity, and so the location 
of the impurity electronic levels inside the band gap of 
the pure crystal. This brief introduction highlights the 
importance of having a good understanding of the lattice 
distortions induced by substitutional impurities in ionic 
crystals. Morgcwer, their experimental measurement is a 
difhcult task,0 □ and so theoretical calculations become 
an ideal complement to the experimental studies. 

From the theoretical side, to obtain an accurate char- 
acterization of the local structure around a defect is a 
delicate matter. First of all, the full translational sym- 
metry of the crystal is lost, and so Bloch's theorem can 
not be directly applied. One way to avoid this problem 
is to duplicate a finite region of the crystal around the 
impurity to recover full translational symmetry and ex- 
ploit the computational convenience jjpf Bloch's theorem. 
This is done in supercell techniques.Ll The supercell size 
has to be chosen sufficiently large so that the interac- 
tion between the local region influenced by the defect 
and its periodical images do not interact, because that 
interaction would be nonphysical. If one forgets about 
Bloch's theorem, one is left with the cluster approach, 
in which the doped crystal is modeled by a finite clus- 
ter centered on the impurity and embedded in a field 
representing the rest of the host lattice. This cluster ap- 
proach is the one chosen in the present study, and has 



been used in the past to study thuei^ometrical and opti- 
cal properties of doped crystalsJju In this approach it 
is important to achieve an accurate description of both 
the active cluster and the environment. Moreover, those 
two descriptions should be consistent with each other. 
Usual deficiencies found in previous works employing the 
cluster approach are the following: (a) an incomplete 
representation of the environment. In the simplest and 
most frequently used approach, the environment is sim- 
ulated by placing point charges on the lattice sites. This 
procedure has to be improved in order to obtain a re- 
alistic desLcpi}tipipjafJ,hepk^^ distortions around the 
impurityBiimemeMEl Model potentials have 
been developed to represent the effects of the environ- 
ment on the active cluster, that include attractive and 
repulsive quantum-mechanical terms aside from the clas- 
sical Madelung term.O (b) An active cluster size that is 
too small. This is usually a problem of the computa- 
tionally expensive ah initio calculations. Only the posi- 
tions of the ions in the first shell araund.|t.hp .irnpurity 
are allowed to relax in most cases. EJoO'EZIljO How- 
ever, geometrical relaxations far beyond the first shell of 
neighbors can be expected for certain impurities, as-sufc 
gested by recent semiempirical simulations of solids.EjCJ 
(c) An abrupt connection between the active cluster and 
the environment. The wavefunctions and positions of 
the ions in the environment are not allowed to relax, so 
that the ions in the surface of the active clustePnfejsl a 
wrong embedding potential. In previous works pHr^ we 
have shown how an unphysically abrupt connection will 
lead to wrong distortions independently of the intrinsic 
accuracy of the electronic structure code, (d) An in- 
complete self-embedding consistency. The cluster model 
must be tested with calculations on the pure crystals in 
order to sufujess systematic errors from the calculated 
distortions .EiJo 

The most expensive ah initio calculations become im- 
practical when the defect induces large-scale lattice dis- 
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tortions. On the other hand, semiempirical calcula- 
tions employing pair potentials are not completely re- 
liable due to the problem of transferability to environ- 
ments different from those in which they were gener- 
ated, la this paper we employ the ab initio Perturbed Ion 
(aiPI)ETEamodel to study the lattice distortions induced 
by Ga+, In+ and T1+ impurities in NaCl, NaBr, KCl and 
KBr. With the use of this method we circunvent all the 
incoveniences listed above: (a) The active cluster is em- 
bedded in a quantum environment represented by the ab 
initio model potentials of Huzinaga et aZHj; (b) The com- 
putational simplicity of the PI model allows for the geo- 
metrical rela xaUo i ii o£^everal coordination shells around 
the impurity (c) The local region around the 
defect in which structural rearrangements are important 
is connected to the frozen crystalline environment by a 
smooth interface of ions whose wavefunctions are allowed 
to selfconsistently relax; (d) Parallel cluster model calcu- 
lations on the pure crystals are performed in order to ob- 
tain trustworthy distortions. By studying several related 
systems we try to identify systematic distortion trends 
that might be useful in later theoretical studies of doped 
crystals sirailar to those here considered. In a previous 
publications we considered the case of ns^ anionic sub- 
stitutional impurities. Now we complete the study of ns^ 
substitutional impurities in alkali halide crystals with a 
consideration of cationic impurities. 

The remainder of this paper is organized as follows: 
In Section O we describe the active clusters which have 



been used to model the doped systems. In Section III we 
present and discuss the results of the calculations, and 
Section IV summarizes the main conclusions. 



II. CLUSTER MODEL 

The ab initio Perturbed Ion model is a particular ap- 
plication of the theory of electronic separability of Huzi- 
naga and coworkerscBE2l to ionic solids in which the basic 
building blocks are reduced to single ions.. The PI model 
was first developed for perfect crystals .lj Its application 
to the study of im purity centers in ionic crystals has been 
described in refs. 15|,|l8|,^ and we refer to those papers 
for a full account of the method. In brief, an active cluster 
containing the impurity is considered, and the Hartree- 
Fock-Roothaan (HFR) equationsEil for each ion in the ac- 
tive cluster are solved in the field of the other ions. The 
Fock operator includes, apart from the usual intra-atomic 
terms, an accurate quantum-mechanical crystal poten- 
tial and a lattice projection operator which accounts for 
the energy contribution due_to the overlap between the 
wave functions of the ionsEj The atomic-like HFR so- 
lutions are used to describe the ions in the active clus- 
ter in an iterative stepwise procedure. The wave func- 
tions of the lattice ions external to the active cluster 
are taken from a PI calculation for the perfect crystal 
and are kept frozen during the embedded-cluster calcu- 



lation. Those wave functions are explicitely considered 
for ions up to a distance d from the center of the active 
cluster such that the quantal contribution from the most 
distant frozen shell to the effective cluster energy is less 
than 10"^ hartree. Ions at distances beyond d contribute 
to the effective energy of the active cluster just through 
the long-range Madelung interaction, so they are repre- 
sented by point charges. At the end of the calculation, 
the ionic wave functions are selfconsistent within the ac- 
tive cluster and consistent with the frozen description of 
the rest of the lattice. The intraatomic Coulomb correla- 
tion, which is neglected at the Hartree- Fock level, is com- 
puted as a correction by usipa-the Coulomb-Hartree-Fock 
(CHF) model of Clementi.BcJ The polarization contri- 
bution is computed by using the p plari zable ion model 
devised by Madden aiid_f;oworkers,L3E3 as explained in 
previous publications Ej'lj In-crystal polarizabilities were 
obtained from Ref. |5^. 

Now we describe the cluster models of increasing com- 
plexity employed to simulate the impurity systems. All of 
them are embedded in accurate quantum environments 
as indicated in the paragraph above. The active clusters 
are split up into two subsets, which we call Ci and C2 
following Refs. Both the positions and wavefunc- 

tions of the ions in the inner Ci subset are allowed to 
relax. The positions of the ions in the outer C2 sub- 
set are held fixed during the optimization process, but 
their wavefunctions are selfconsistently optimized. Thus, 
the C2 subset provides a smooth interface connecting 
the Ci region, where distortions are important, to the 
frozen crystalline environment. In practice, C2 contains 
all those ions which are first neighbors of the ions in Ci 
and are not already contained in Ci. In Table I we show 
the lattice ions included in the Ci subset for the different 
cluster models. It has to be understood that the impu- 
rity central ion is included in all cluster models, and that 
each cluster model is formed by adding the ions shown 
in Table I to the preceding cluster model. In each case, 
the geometrical relaxation around the impurity has been 
performed by allowing for the independent breathing dis- 
placements of each shell of ions, and minimizing the total 
energy with respect to those displacements until the ef- 
fective energies aije-, converged up to 1 meV. A downhill 
simplex algorithmEJ was used. For the description of the 
ions we have uspi large STO basis setSj-all taken from 
Clementi-RoettiEil and McLean-McLeanN tables. The 
calculations have been perfacmed by employing the ex- 
perimental lattice constanta£j to describe the geometri- 
cally frozen part of the crystals. 

The self-embedding consistency of the method has 
been checked for all the different cluster models and 
found to be of .a quality similar to that found in previ- 
ous papers.EIlS By this we mean that if the pure crystal 
is represented by one of the cluster models enumerated 
above (that is, if the central impurity is replaced by the 
alkali cation corresponding to the pure crystal), the re- 
sults of the cluster model calculations closely agree with 
those from a PI calculation for the pure crystal, where all 
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cations (or anions) are equivalent by translational sym- 
metry. Nevertheless, the self-embedding consistency is 
never complete. In order to supress systematic errors 
from the distortions calculated with the cluster method, 
the radial displacements of each shell have been calcu- 
lated using the following formula: 

ARi^ R,{Imp+ : AX)~ Ri{A+ : AX), (1) 

where Ri (i=l, 2, 3, 4, 7, 8) refer to the radii of the first, 
second, third, fourth, seventh and eigth shells around 
the impurity in the AX crystal, A = Na, K, X = CI, Br, 
and Imp+ — Ga"*", In"'", Tl"*". Thus both systems (pure 
and doped crystals) are treated in eq. (1) on equal foot 
so that any inaccuracy in the calculated distortions is 
completely due to the electronic structure method. 

III. RESULTS AND DISCUSSION 
A. Influence of the active cluster size 

We show in Table II the lattice distortions calculated 
with the different cluster models described in the previ- 
ous section, for the NaChTl"*" system. As the distortion 
trends for different impurity systems will be discussed 
in the next subsection, here we just point out two gen- 
eral results; (a) the distortion in the first coordination 
shell, Ai?i, depends upon an explicit consideration of the 
distortion in the fourth coordination shell, Ai?4, that is 
the distortion along the (100) crystallographic direction 
propagates at least until next nearest neighbor positions. 
This is clearly seen in Table II by comparing the value of 
Ai?i in cluster model C with the value of Ai?i in cluster 
models D and E. When the positions of the ions in the 
fourth coordination shell are held fixed, the Ai?i expan- 
sion can not attain its optimal value; (b) the distortions 
along (100) and (110) crystallographic directions arc cou- 
pled. This means that the calculated distortions in the 
first and fourth coordination shells depend upon an ex- 
plicit inclusion of the second and eigth coordination shells 
in the active clusters. This is seen in the value of ARi 
when passing from model A to model B and in Ai?i and 
Ai?4 when passing from model E to model F. These con- 
clusions are independent of the specific system studied. 
For example, we show in Fig. 1 the evolution of Ai?i, as 
a function of the active cluster complexity, for the three 
cationic impurities in NaCl. We may appreciate that the 
three lines are neatly parellel to each other, supporting 
the validity of our conclusions. Thus, the errors in Ai?i, 
due only to a small active cluster size (that is not directly 
associated with the accuracy of the electronic structure 
code), can be roughly estimated to be 1%. 

A magnitude which is very difficult to estimate theo- 
retically is the defect formation energy. At low pressure 
and temperature conditions, the formation of the impu- 
rity centers should be discussed in terms of-ithe internal 
energy difference for the exchange reactionlHJ 



{A+ : AX)s + Imp+ ^ {Imp+ : AX)s + A+ , (2) 

where the s and g subindexes refer to solid and gas 
phases, respectively. The Madelung energy term, which 
is the most important contribution to cohesion in ionic 
crystals, is minimized when the ions are in their perfect 
crystallographic positions, so the lattice distortions de- 
scribed above will tend to destabilize this term. The de- 
fect will be stable only if the electronic contributions can 
compensate for that Madelung energy loss. The detailed 
energy balance contains then nonnegligible contributions 
from energy terms, like the dispersion interactions, that 
arc not so important for the determination of structural 
properties. The dispersion interactions, as well as any 
relativistiC|-energy terms, are not included in the aiPI 
formalism,E3 so wc do not expect the formation energies 
calculated from it to be wholly reliable. Nevertheless, 
the contribution to the defect formation energy coming 
from elastic relaxation of the lattice is meaningful, and 
we show its evolution as a function of cluster complex- 
ity in Fig. 2, for the case of Ga"*", In"*" and T1+ im- 
purities in NaCl. Specifically, we show the defect forma- 
tion energy differences (stabilization energies), calculated 
by substracting the formation energies calculated with a 
given cluster model and its predecesor. As a reference 
state for the cluster model A we took a smaller active 
cluster containing just 7 ions, namely the impurity and 
its 6 nearest neighbor Cl~ anions, which were held fixed 
at their lattice positions. This cluster can be considered 
the limit case of those shown in the previous section, 
where the Ci subset is formed by the impurity cation 
and the C2 subset (the interphase) is formed just by the 
six first neighbors of the impurity. The figure seems to 
suggest that the first shell distortions sligtly destabilize 
the defects, while second shell distortions largely stabi- 
lize them. This apparently paradoxical result is just a 
consequence of our stepped construction of cluster mod- 
els. If we include the second shell distortions alone, the 
net effect is also a destabilization. This simply shows 
that the distortions in the first and second shells, where 
short-range overlap interactions with the impurity are 
important, should be included together in the calcula- 
tion. Structural relaxation of those two shells provide 
the largest contribution to elastic relaxation. We also 
appreciate nonnegligible stabilization energies when the 
third, fourth, and eigth coordination shells are allowed 
to relax, giving an idea of the difficulty of obtaining well 
converged theoretical formation energies. 

B. Distortion trends 

In this section we show the calculated distortions in- 
duced by substitutional Ga"*", In"*" and Tl"*" impurities in 
NaCl, NaBr, KCl and KBr crystals. For all these calcu- 
lations we employ the cluster model F, which is the most 
complete one from those considered in previous sections. 
These distortions, collected in Table III, are the main 
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quantitative result from our study. From that table we 
can extract interesting distortion trends, that we describe 
in the following: 

First shell distortion. This shell undergoes an expan- 
sion in all cases, with the only exception ofj-the Ga^ im- 
purity in K salts. In a previous publicationtJ we showed 
how these distortions can be rationalized in terms of 
simple geometric arguments involving ion sizes. As a 
rough measure of ion size we employed the quantity r 
= (< >y^^, where the expectation value is taken 
over the outermost orbital of the ion, and is calculated 
from the crystal-consistent ionic wave functions obtained 
through the aiPI calculations. Then the first shell dis- 
tortions were shown to be related to the differences 6 
=r(Imp+:AX)-r(A+:AX), where r(Imp+:AX) is the ra- 
dius of the Ga"*", In+ or T1+ impurities in the AX crys- 
tal, and r(A+:AX) is the radius of the alkali cation in 
the pure AX crystal. A similar analysis here shows that 
6 is negative only for KCl:Ga+ and KBr:Ga+, giving an 
explanation to the inwards relaxation of the first shell in 
those two cases. More in general, for a fixed crystal 5 
increases in the series Ga"*" — > In"*" T1+, as should be 
expected. The ARi distortions for any given crystal fol- 
low the same trend. For a given impurity and fixed halide 
anion, the first shell expansion is larger the smaller the 
size of the alkali cation. This trend suggests what kind of 
distortions can be expected in Li or Rb host lattices. For 
example, one can expect a small inwards relaxation for 
RbCl:In+, similar to that found for KCl:Ga+, and more 
important distortions for Li salts, that perhaps will ask 
for larger active cluster sizes. For those systems where 
an expansion is observed, if we fixed the impurity cation 
and the alkali cation, the expansion is larger the larger 
the size of the anion. Thus one can expect smaller expan- 
sions for NaF and KF salts and larger ones for Nal and 
KI salts. For those cases where a contraction is observed, 
the contraction is larger the smaller the size of the anion. 

Second shell distortion. The displacement of the sec- 
ond shell is always an expansion. For any given host 
lattice, that expansion is larger the larger the impurity 
size, because of the increase in cation-cation overlap re- 
pulsion. If the impurity and the halide anion are fixed, 
the expansion increases with decreasing alkali size. This 
trend suggests again larger Ai?2 expansions in Li salts 
and smaller ones in Rb salts. If the impurity and the 
alkali cation are fixed, the expansion is larger the smaller 
the anion size, because cation-cation overlap decreases 
with increasing anion size. Then, larger AR2 expansions 
are expected in NaF and KF salts and smaller ones for 
Nal and KI salts. 

Third shell distortion. This shell experiences a small 
contraction in all cases. For a given crystal, the con- 
traction is larger the larger the size of the impurity. If 
we fix the impurity and the halide anion, the contrac- 
tion increases with decreasing alkali size. Then, larger 
contractions should be expected in Li salts, and smaller 
ones in Rb salts. If the impurity and the alkali cation 
are fixed, the contraction is larger the smaller the size of 



the halide anion, so larger contractions are expected for 
fluoride salts. As the direct overlap interactions between 
the impurity and the ions in the third coordination shell 
are negligibly small, this contraction has to be consid- 
ered an indirect doping effect, mediated by the distor- 
tions induced on the first and second shells. The small 
contraction of the third shell optimizes the Madelung en- 
ergy around the impurity and also serves to pack more 
efficiently the ions in response to the outward motion of 
the ions in the first and second shells. 

Fourth and seventh shell distortions. The fourth shell 
experiences an expansion, except in the case of KCl:Ga+. 
The AR4 displacements proceed along the same (100) 
crystallographic direction as the Ai?i displacements, and 
thus are clearly induced by the first shell distortions, 
as indicated in the previous subsection. The distortion 
trends of this shell are then found to be the same as those 
of the first shell. The displacements of those ions in the 
seventh coordination shell are always negligibly small, in- 
dicating that the distortions propagate along the (100) 
direction just until next-nearest neighbor positions. 

Eigth shell distortion. Similarly to the previous para- 
graph, the eigth shell undergoes an expansion induced by 
the expansion of the second shell, as both displacements 
proceed along the (110) crystallographic directions. Thus 
the distortion trends are the same as those found for the 
second shell. 



IV. SUMMARY 

A study of the local lattice distortions induced by sub- 
stitutional Ga+, In+ and T1+ impurities in NaCl, NaBr, 
KCl and KBr crystals has been reported. For that pur- 
pose we have considered active clusters of increasing com- 
plexity, with a number of ions varying between 33 and 
177 ions, embedded in accurate quantum environments 
representing the rest of the crystal. The local distortions 
obtained extend beyond the first shell of neighbors in all 
cases. Thus, the assumptions frequently employed in im- 
purity calculations, which consider the active space as 
formed by the central impurity plus its first coordina- 
tion shell only, should be taken with some care. More- 
over, we have found that the distortions along the (100) 
and (110) crystallographic directions are not independent 
from each other, and also that the first shell distortion 
is not converged if the positions of the ions in the fourth 
coordination shell around the impurity are not allowed 
to relax. Once we have obtained a reliable cluster model, 
distortion trends have been identified and discussed. In 
those cases where the size of the impurity is larger than 
that of the alkali cation, the impurity induces an expan- 
sion propagating along (100) and (110) crystallographic 
directions until next-nearest neighbor positions. There 
is also an indirect contraction along the (111) crystallo- 
graphic direction, that probably affects only to the near- 
est neighbor ions in that direction, and which serves to 
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partly compensate the Madclung field around the impu- 
rity and to pack more efficiently the ions in response to 
the outward motion of the first and second shells. In 
those cases where the size of the impurity is smaller than 
that of the alkali cation of the host lattice, the only sig- 
nificant change in these trends is the inwards relaxation 
along the (100) direction. The dependence of the several 
shell distortions upon a change of alkali cation, halide 
anion, and impurity cation have been described. From 
this description one can advance expected distortions for 
closely related systems. 
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Captions of Tables 

Table I. Description of the lattice ions included in the 
Ci subset for the different cluster models. The degener- 
acy (number of symmetry equivalent ions) of each coor- 
dination shell is shown, together with the charge of the 
ions forming that shell, and the total number of atoms 
in the active cluster (ions in Ci and C2 subsets). 

Table II. Lattice distortions (in A) and percentage 
distortion values obtained with the different cluster mod- 
els described in the text. 

Table III. Lattice distortions (in A) of several coor- 
dination shells around Ga"'", In+ and Tl"*" impurities in 
four alkali halide crystals, employing cluster model F as 
described in the text. 

Captions of Figures 

Figure 1. First-shell distortions (in A) induced by 
Ga+, In+ and T1+ impurities in NaCl, as a function of 
the complexity of the cluster model. Note the similar 
behavior of the three curves. 

Figure 2. Defect formation energy stabilization (in 
eV) as a function of cluster complexity. 
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